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Abstract. The Galactic population of close white dwarf binaries is expected to provide the 
largest number of gravitational wave sources for low frequency detectors such as the Laser 
Interferometer Space Antenna (LISA). Current data analysis techniques have demonstrated the 
capability of resolving on the order of 10* white dwarf binaries from a 2 year observation. 
Oh' Resolved binaries are either at high frequencies or large amplitudes. Such systems are more 
likely to be high-mass binaries, a subset of which will be progenitors of SNe la in the double 
degenerate scenario. We report on results of a study of the properties of resolved binaries using 
' a population synthesis model of the Galactic white dwarf binaries and a LISA data analysis 
^Nj , algorithm using Mock LISA Data Challenge tools. 
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1^ Space-based gravitational wave detectors, such as the Laser Interferometer Space An- 

tenna (LISA) (Danzmann 1996, Hughes 2006), will be sensitive to gravitational radiation 
in the millihertz range. The Galactic population of double white dwarf (DWD) binaries 
C3 , will be the dominant source in this frequency range. With some 30 million systems 
thought to be producing signals with frequencies above 0.1 mHz, these binaries will pro- 
I , duce a confusion-limited signal below 3 mHz, with multiple sources in each resolvable 
J> ' frequency bin. Consequently, most individually resolvable signals will come from DWDs 
OO ] that are of high frequency or high mass, or both. The typical double degenerate progeni- 
tor of a type la supernova consists of two carbon-oxygen white dwarfs in an ultra-compact 
binary. Many of these are in the class of individually resolvable binaries, and so LISA will 
be a double degenerate finder that can identify many of the progenitor systems within the 
Milky Way. In order to estimate any selection biases within the binaries observed with 
LISA, we have generated a two-year LISA data stream using the DWD population from 
Ruiter et al. (2010). We have then applied a detection and parameter estimation algo- 
rithm developed by us (Stroeer & Benacquista 2011) in order to determine the properties 
of the individually resolved binaries. 
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1. Introduction 



2. Analysis 

We model the population of DWDs for the bulge and the disk of the Galaxy using the 
StarTrack population synthesis code of Belczynski (Belczynski et al. 2002, 2008). We 
assume a constant star formation rate over the past 10 Gyr for the disk and a constant 
star formation over the first Gyr for the bulge. This results in a population of nearly 30 
million DWDs within the LISA sensitivity band. The details of the population synthesis 
and distribution can be found in Ruiter et a/. (2010). Using this model, we compute 
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the LISA response to the entire population using the standard LISA configuration of 5 
million km armlengths and a two-year observation time. We then compute the integrated 
signal-to-noise ratio (p) for every DWD in the population using: 



where h is the Fourier transform of the expected gravitational wave strain due to a single 
DWD, and Sn{f) is the one sided instrumental noise spectral density. We select only 
those DWDs with p ^ 5 to do a full search using our data analysis pipeline. This first cut 
gives about 40 thousand binaries. The data analysis pipeline identifies about 31 thousand 
systems. These are the DWDs we describe as resolved. We define a Type la progenitor 
as any DWD with a total mass above 1.44 Mq. In the underlying Galactic population, 
there are 405,309 SNe la progenitors out of 27,835,248 DWDs. This ratio drops to 653 
out of 39,468 after the first cut, and ends at 503 out of 31,169 resolved systems. 

The number of resolved systems are heavily weighted toward systems with gravitational 
wave frequencies above 1 mHz. This can be explained by the fact that there are more 
than several hundred binaries in each resolvable frequency below 1 mHz, and so systems 
with lower frequencies are more likely to be lost in the confusion of signals from all the 
binaries in the same frequency bin. Above this frequency, nearly all systems are recovered, 
as can be seen in Figure [TJ Thus, the relative number of SNe la progenitors compared to 
the total population can be assumed to be complete at higher frequencies. 

When we look at the distribution of distances to the resolved binaries compared to the 
underlying Galactic population, we find that there is no bias towards nearby systems, 
as can be seen in Figure [2l This can be understood by realizing that any binary with 
an orbital frequency above the confusion-limited cut-off of about 1 mHz will be bright 
enough in gravitational waves to be visible across the Galaxy. However, as \h\ oc d~^, we 
do not detect distant systems with d > 30 kpc. 

Parameter estimation in LISA for DWDs returns the amplitude {\h\), the sky location 
and orientation, and (if measurable) the "chirp" or frequency shift due to the emission 
of gravitational radiation as well as mass transfer and tidal interactions if applicable. In 
order to identify which resolved binaries are SNe la progenitors, we need to determine 
the total mass. The "chirp mass", is given by Mchirp = fM^^^M^^^ where /i is the reduced 
mass and M is the total mass. If the chirp (/) can be measured, then Mchirp can be 
found through Mchirp fx f^^^ f^^^^^ , where the proportionality only depends on physical 
constants. Once Af chirp is measured and the system is "clean" so that only gravitational 
radiation contributes to the chirp, then the distance to the binary can be obtained (Schutz 
1997), and so the spatial distribution of these systems throughout the Galaxy will be 
known. The chirp masses of the resolved binaries are shown in Figure [31 and it can be 
seen that all resolved systems with Mchirp ^ 0.62 Mq are SNe la progenitors. The spatial 
distribution of the resolved binaries is shown in Figure |4] and indicates that they are 
found throughout the Galaxy. 

3. Conclusions 

We have found that LISA will be able to resolve less than 0.1% of the Galactic DWDs 
within its sensitivity band of / ^ 0.1 niHz. In the underlying population of our synthetic 
Galaxy, 1.4% of these DWDs are potential SNe la progenitors. The resolved population 
consists of nearly every binary with a gravitational wave frequency above 3 mHz, and 
that this population is observable throughout the Galaxy. Type la progenitors make 
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Figure 1. Number in frequency space for tlie underlying population of double white dwarfs (top 
panel), for systems with a signal-to- noise ratio greater than 5 (middle panel), and for systems 
that were resolved using parameter estimation (bottom panel). 
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Figure 2. Histograms of distances to double white dwarfs. The upper panel shows all binaries 
while the lower panel is only potential Type la progenitors. The underlying population is in black, 
the systems with a signal-to-noise ratio greater than 5 are in red, and the resolved systems are 
shown in green. 

up about 1.4% of the resolved systems, and so there are no appreciable biases in this 
population. Thus, LISA can be considered to be an effective tool for identifying the 
Galactic population of DWDs and the SNe la progenitors within it. 
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Figure 3. Distribution of chirp masses for the underlying population of double white dwarfs 
and for the population of resolved binaries with measurable chirp (defined as /Tobs > !)■ The 
underlying population is shown in the solid line. The population of all resolved binaries is the 
dashed line, and the population with measurable chirp is indicated by the filled area. Type la 
progenitors are shown in red. 




Figure 4. Spatial distribution of resolved binaries viewed from above the Galaxy in galactocen- 
tric coordinates. All binaries are shown in black while potential Type la progenitors are shown 
in red. The approximate location of the sun is shown in cyan. 



Danzmann, K. 1996, CQG, 13, 247 

Hughes, S. 2006, in S.M. Merkowitz & J.C. Livas (eds). Laser Interferometer Space Antenna: 

6th Int. LISA Symp. (Melville, NY: AIP) p. 13 
Ruiter, A. J., Belczynski, K., Benacquista, M., Larson, S.L., & Williams, G. 2010, ApJ, 717, 1006 
Schutz, B. 1997, in Proceedings of the Alpbach Summer School 1997, 'Fundamental Physics in 

Space' (ESA-SP-420) p. 229 
Stroeer, A. & Benacquista, M. 2011, in preparation 



